Transcriptomic analysis of the human habenula in schizophrenia.

Importance: Habenula (Hb) pathophysiology is involved in many neuropsychiatric disorders, including schizophrenia. Deep brain stimulation and pharmacological targeting of the Hb are emerging as promising therapeutic treatments. However, little is known about the cell type-specific transcriptomic organization of the human Hb or how it is altered in schizophrenia. Objective: To define the molecular neuroanatomy of the human habenula and identify transcriptomic changes in individuals with schizophrenia compared to neurotypical controls. Design Setting and Participants: This study utilized Hb-enriched postmortem human brain tissue. Single nucleus RNA-sequencing (snRNA-seq) and single molecule fluorescent in situ hybridization (smFISH) experiments were conducted to identify molecularly defined Hb cell types and map their spatial location (n=3-7 donors). Bulk RNA-sequencing and cell type deconvolution were used to investigate transcriptomic changes in Hb-enriched tissue from 35 individuals with schizophrenia and 33 neurotypical controls. Gene expression changes associated with schizophrenia in the Hb were compared to those previously identified in the dorsolateral prefrontal cortex (DLPFC), hippocampus, and caudate. Main Outcomes and Measures: Semi-supervised snRNA-seq cell type clustering. Transcript visualization and quantification of smFISH probes. Bulk RNA-seq cell type deconvolution using reference snRNA-seq data. Schizophrenia-associated gene differential expression analysis adjusting for Hb and thalamus fractions, RNA degradation-associated quality surrogate variables, and other covariates. Cross-brain region schizophrenia-associated gene expression comparison. Results: snRNA-seq identified 17 cell type clusters across 16,437 nuclei, including 3 medial and 7 lateral Hb populations. Cell types were conserved with those identified in a rodent model. smFISH for cell type marker genes validated snRNA-seq Hb cell types and depicted the spatial organization of subpopulations. Bulk RNA-seq analyses yielded 45 schizophrenia-associated differentially expressed genes (FDR < 0.05), with 32 (71%) unique to Hb-enriched tissue. Conclusions: These results identify topographically organized cell types with distinct molecular signatures in the human Hb. They further demonstrate unique transcriptomic changes in the epithalamus associated with schizophrenia, thereby providing molecular insights into the role of Hb in neuropsychiatric disorders.

Single-Nucleus Transcriptome Profiling of Dorsolateral Prefrontal Cortex: Mechanistic Roles for Neuronal Gene Expression, Including the 17q21.31 Locus, in PTSD Stress Response.

OBJECTIVE: Multidisciplinary studies of posttraumatic stress disorder (PTSD) and major depressive disorder (MDD) implicate the dorsolateral prefrontal cortex (DLPFC) in disease risk and pathophysiology. Postmortem brain studies have relied on bulk-tissue RNA sequencing (RNA-seq), but single-cell RNA-seq is needed to dissect cell-type-specific mechanisms. The authors conducted the first single-nucleus RNA-seq postmortem brain study in PTSD to elucidate disease transcriptomic pathology with cell-type-specific resolution. METHOD: Profiling of 32 DLPFC samples from 11 individuals with PTSD, 10 with MDD, and 11 control subjects was conducted (∼415K nuclei; >13K cells per sample). A replication sample included 15 DLPFC samples (∼160K nuclei; >11K cells per sample). RESULTS: Differential gene expression analyses identified significant single-nucleus RNA-seq differentially expressed genes (snDEGs) in excitatory (EX) and inhibitory (IN) neurons and astrocytes, but not in other cell types or bulk tissue. MDD samples had more false discovery rate-corrected significant snDEGs, and PTSD samples had a greater replication rate. In EX and IN neurons, biological pathways that were differentially enriched in PTSD compared with MDD included glucocorticoid signaling. Furthermore, glucocorticoid signaling in induced pluripotent stem cell (iPSC)-derived cortical neurons demonstrated greater relevance in PTSD and opposite direction of regulation compared with MDD, especially in EX neurons. Many snDEGs were from the 17q21.31 locus and are particularly interesting given causal roles in disease pathogenesis and DLPFC-based neuroimaging (PTSD: ARL17B, LINC02210-CRHR1, and LRRC37A2; MDD: LRRC37A and LRP4), while others were regulated by glucocorticoids in iPSC-derived neurons (PTSD: SLC16A6, TAF1C; MDD: CDH3). CONCLUSIONS: The study findings point to cell-type-specific mechanisms of brain stress response in PTSD and MDD, highlighting the importance of examining cell-type-specific gene expression and indicating promising novel biomarkers and therapeutic targets.

Orthogonal approaches required to measure proteasome composition and activity in mammalian brain tissue.

Proteasomes are large macromolecular complexes with multiple distinct catalytic activities that are each vital to human brain health and disease. Despite their importance, standardized approaches to investigate proteasomes have not been universally adapted. Here, we describe pitfalls and define straightforward orthogonal biochemical approaches essential to measure and understand changes in proteasome composition and activity in the mammalian central nervous system. Through our experimentation in the mammalian brain, we determined an abundance of catalytically active proteasomes exist with and without a 19S cap(s), the regulatory particle essential for ubiquitin-dependent degradation. Moreover, we learned that in-cell measurements using activity-based probes (ABPs) are more sensitive in determining the available activity of the 20S proteasome without the 19S cap and in measuring individual catalytic subunit activities of each ß subunit within all neuronal proteasomes. Subsequently, applying these tools to human brain samples, we were surprised to find that post-mortem tissue retained little to no 19S-capped proteasome, regardless of age, sex, or disease state. In comparing brain tissues (parahippocampal gyrus) from patients with Alzheimer's disease (AD) and unaffected individuals, the available 20S proteasome activity was significantly elevated in severe cases of AD, an observation not previously noted. Taken together, our study establishes standardized approaches for the comprehensive investigation of proteasomes in mammalian brain tissue, and we reveal new insight into brain proteasome biology.

SUFI: an automated approach to spectral unmixing of fluorescent multiplex images captured in mouse and post-mortem human brain tissues.

BACKGROUND: Multispectral fluorescence imaging coupled with linear unmixing is a form of image data collection and analysis that allows for measuring multiple molecular signals in a single biological sample. Multiple fluorescent dyes, each measuring a unique molecule, are simultaneously measured and subsequently "unmixed" to provide a read-out for each molecular signal. This strategy allows for measuring highly multiplexed signals in a single data capture session, such as multiple proteins or RNAs in tissue slices or cultured cells, but can often result in mixed signals and bleed-through problems across dyes. Existing spectral unmixing algorithms are not optimized for challenging biological specimens such as post-mortem human brain tissue, and often require manual intervention to extract spectral signatures. We therefore developed an intuitive, automated, and flexible package called SUFI: spectral unmixing of fluorescent images. RESULTS: This package unmixes multispectral fluorescence images by automating the extraction of spectral signatures using vertex component analysis, and then performs one of three unmixing algorithms derived from remote sensing. We evaluate these remote sensing algorithms' performances on four unique biological datasets and compare the results to unmixing results obtained using ZEN Black software (Zeiss). We lastly integrate our unmixing pipeline into the computational tool dotdotdot, which is used to quantify individual RNA transcripts at single cell resolution in intact tissues and perform differential expression analysis, and thereby provide an end-to-end solution for multispectral fluorescence image analysis and quantification. CONCLUSIONS: In summary, we provide a robust, automated pipeline to assist biologists with improved spectral unmixing of multispectral fluorescence images.

Integrated DNA methylation and gene expression profiling across multiple brain regions implicate novel genes in Alzheimer's disease.

Late-onset Alzheimer's disease (AD) is a complex age-related neurodegenerative disorder that likely involves epigenetic factors. To better understand the epigenetic state associated with AD, we surveyed 420,852 DNA methylation (DNAm) sites from neurotypical controls (N = 49) and late-onset AD patients (N = 24) across four brain regions (hippocampus, entorhinal cortex, dorsolateral prefrontal cortex and cerebellum). We identified 858 sites with robust differential methylation collectively annotated to 772 possible genes (FDR < 5%, within 10 kb). These sites were overrepresented in AD genetic risk loci (p = 0.00655) and were enriched for changes during normal aging (p < 2.2 × 10-16), and nearby genes were enriched for processes related to cell-adhesion, immunity, and calcium homeostasis (FDR < 5%). To functionally validate these associations, we generated and analyzed corresponding transcriptome data to prioritize 130 genes within 10 kb of the differentially methylated sites. These 130 genes were differentially expressed between AD cases and controls and their expression was associated with nearby DNAm (p < 0.05). This integrated analysis implicates novel genes in Alzheimer's disease, such as ANKRD30B. These results highlight DNAm differences in Alzheimer's disease that have gene expression correlates, further implicating DNAm as an epigenetic mechanism underlying pathological molecular changes associated with AD. Furthermore, our framework illustrates the value of integrating epigenetic and transcriptomic data for understanding complex disease.

Genetic risk mechanisms of posttraumatic stress disorder in the human brain.

Posttraumatic stress disorder (PTSD) follows exposure to a traumatic event in susceptible individuals. Recently, genome-wide association studies have identified a number of genetic sequence variants that are associated with the risk of developing PTSD. To follow up on identifying the molecular mechanisms of these risk variants, we performed genotype to RNA sequencing-derived quantitative expression (whole gene, exon, and exon junction levels) analysis in the dorsolateral prefrontal cortex (DLPFC) of normal postmortem human brains. We further investigated genotype-gene expression associations within the amygdala in a smaller independent RNA sequencing (Genotype-Tissue Expression [GTEx]) dataset. Our DLPFC analyses identified significant expression quantitative trait loci (eQTL) associations for a "candidate" PTSD risk SNP rs363276 and the expression of two genes: SLC18A2 and PDZD8, where the PTSD risk/minor allele T was associated with significantly lower levels of gene expression for both genes, in the DLPFC. These eQTL associations were independently confirmed in the amygdala from the GTEx database. Rs363276 "T" carriers also showed significantly increased activity in the amygdala during an emotional face-matching task in healthy volunteers. Taken together, our preliminary findings in normal human brains represent a tractable approach to identify mechanisms by which genetic variants potentially increase an individual's risk for developing PTSD. © 2016 Wiley Periodicals, Inc.

Implementation and clinical characteristics of a posttraumatic stress disorder brain collection.

A postmortem human brain collection to study posttraumatic stress disorder (PTSD) is critical for uncovering the molecular mechanisms that contribute to this psychiatric disorder. We describe here the PTSD brain collection at the Lieber Institute for Brain Development in Baltimore, Maryland, consisting of postmortem brain donations acquired between 2012 and 2017. Thus far, 87 brains from individuals meeting DSM-5 criteria for PTSD were collected after consent was obtained from legal next-of-kin, and subsequently clinically characterized for molecular studies. PTSD brain donors had high rates of comorbid diagnoses, including depression (62.1%), substance abuse (74.7%), drug-related death (69.0%), and suicide completion (17.2%). PTSD cases were subdivided into two categories: combat-related PTSD (n = 24) and noncombat/domestic PTSD (n = 63). The major differences between the combat-related and domestic PTSD cohorts were sex, drug-related death, and the prevalence of bipolar disorder (BPD) comorbidity. The combat-related group was entirely male, with only one BPD subject (4.2%), and had significantly fewer drug-related deaths (45.8%) in contrast to the domestic group (31.8% male, 36.5% bipolar, and 77.8% drug-related deaths). Medical examiners' offices, particularly in areas with higher military populations, are an excellent source for PTSD brain donations of both combat-related and domestic PTSD.